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The evolution of the entanglement between oscillators that interact with the same environment 
displays highly non-trivial behavior in the long time regime. When the oscillators only interact 
through the environment, three dynamical phases were identified (J. P. Paz and A. Roncaglia, Phys. 
Rev. Lett. 100 (2008)) and a simple phase diagram characterizing them was presented. Here we 
generalize those results to the cases where the oscillators are directly coupled and we show how a 
degree of mixidness can affect the final entanglement. In both cases, entanglement dynamics is fully 
characterized by three phases (SD: sudden death, NSD: no-sudden death and SDR: sudden death 
and revivals) which cover a phase diagram that is a simple variant of the previously introduced one. 
We present results when the oscillators are coupled to the environment through their position and 
also for the case where the coupling is symmetric in position and momentum (as obtained in the 
RWA). As a bonus, in the last case we present a very simple derivation of an exact master equation 
valid for arbitrary temperatures of the environment. 



I. INTRODUCTION 

In recent years it became clear that the study of the 
evolution of entanglement for open quantum systems is 
an important issue not only for fundamental reasons but 
also for practical ones, as entanglement is an essential 
resource for quantum information processing ^ . Entan- 
glement for systems of continuous variable is at the heart 
of the EPR argument [2\ and was discussed in the con- 
text of recent experimental demonstrations of quantum 
teleportation [3] and cryptographic protocols [4 . The 
decoherence induced by the interaction with the envi- 
ronment is an important issue to consider in this con- 
text. In general, decoherence produces dis-entanglement, 
which may occur in a finite time. This phenomenon 
[5J [SJ [7] is known as "sudden death" of entanglement 
(SD). But the fate of entanglement for a quantum open 
system is not at all evident and many authors studied it 
obtaining rather surprising results for systems of qubits 
[5J ini Uni EJ HH HS] and for continuous variable systems 

[iiiTiiiiniiniiiHjiiniiiniiiiiisiiiia. 

More recently, [23J|2S], we provided a unified picture 
that enabled us to understand the various qualitatively 
different types of evolutions of the entanglement in non- 
Mar kovian environments. In fact, we showed that the 
asymptotic dynamics of entanglement can be described 
by three possible phases: SD (sudden death), SDR (sud- 
den death and revival) and NSD (no sudden death) . The 
fate of entanglement can be described using a simple 
phase diagram whose boundaries can even be analytical 
studied in some simple cases. In the above mentioned pa- 
pers we presented the phase diagram under some simple 
assumptions: In particular, we assumed that the oscilla- 
tors did not interact directly (only through the environ- 
ment). Here, we briefly review the results of [24i ,25. and 
we generalize them in two simple ways: we consider more 
general initial conditions and we also consider the case 
when the two oscillators directly interact between them. 

The paper is organized as follows. In Section |ll] we re- 
view the basic notions of entanglement for two harmonic 



oscillators prepared in a Gaussian state. Here, we also 
describe the way in which the evolution of the entan- 
glement can be studied and present a simple quantum 
optical analogy that enables us to understand the qual- 
itative behavior displayed by the entanglement for long 
times. In Section |III| we discuss the models for system- 
environment coupling and present the phase diagram for 
the standard case briefly reviewing the results of [Mll25j . 



In Section IV we present the phase diagram for the case 
where there is a degree of impurity in the state of the vir- 
tual oscillators and direct interactions between the real 
oscillators. Finally, we summarize in Section [V] The Ap- 
pendix[X]contains a simple derivation of the exact master 
equation for an oscillator coupled with a bosonic envi- 
ronment through an interaction term which is symmetric 
under position and momentum interchange (similar to 
what is obtained under the usual RW approximation). 



II. ENTANGLEMENT BETWEEN TWO 
OSCILLATORS 

We will consider a system of two identical quantum 
harmonic oscillators (with coordinates xi and 2:2). The 
interaction with the environment will be discussed in the 
next Section. Here we will discuss how can the entangle- 
ment between such oscillators be quantifled and studied. 
We will restrict to consider Gaussian states (that will 
remain Gaussian under evolution according to the mod- 
els described below) and factorized initial conditions be- 
tween system and environment. For such class of states 
we will be able to obtain simple analytical results. 

Entanglement for Gaussian states of two bosonic 
modes is entirely determined by the covariance matrix 
deflned as 



V,,{t) = {{r„r,})/2-{r.,){r,), 



(1) 



where i,j — 1, ... ,4 and r — (xi,pi, X2,P2)- In fact, a 
good measure of entanglement is the logarithmic nega- 
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tivity i;^- [IS [271 [55] defined as: 

Ej^ = max{0, - ln(2z^min)}, 



(2) 



where ^'min is the smallest symplectic eigenvalue of the 
partially transposed covariance matrix. Therefore, the 
entanglement between the two oscillators is entirely de- 
termined by the second moments contained in Vij. Thus, 
the strategy to study the evolution of the entanglement is 
obvious: We should find out the evolution of the elements 
of such matrix. In the following section we will show how 
to do this in a simple, but realistic, model. But it is use- 
ful to advance here some of the most important results 
as they turn out to be rather independent of the details 
but only on the following important assumption: We will 
consider situations in which the system-environment in- 
teraction only involves a bilinear coupling in the collec- 
tive coordinate x+ = xi+X2 (i-e., the relative coordinate 
X- = xi — X2 will be effectively decoupled from the en- 
vironment). In such cases, the asymptotic state will be 
such that the a;+ oscillator will reach an equilibrium state 
characterized by the dispersions Aa:+ and ■ Such dis- 
persions depend upon parameters of the model (spectral 
densities, initial temperature, etc) and such dependence 
will be discussed below. For the moment we only need 
to assume the existence of an equilibrium state for xj^. 
Also, as equilibrium approaches, correlations between 
and x_ vanish. On the other hand, the second moments 
of X- corresponds to a free oscillator with a certain fre- 
quency UJ-. 

The above simple observations are almost all we need 
to fully analyze the evolution of the entanglement be- 
tween the two oscillators. Thus, in the long time limit, 
the covariance matrix has a simple block-diagonal from 
in terms of the variances of the x± oscillators. From this, 
it is possible to obtain covariances for the original xi^2 
oscillators and to find the smallest symplectic eigenvalue 
of the partially transposed version of such matrix. The 
result for the logarithmic negativity is: 



SaaW ^max{0,£;(i)}, 
where the function E{t) is defined as 

E{t) = EM + ^EMG{t). 



(3) 



(4) 



Here G{t) is an oscillatory function with period tt/lo^ 
that takes values in the interval {—1,+!}. The mean 
value i?7v and the amplitude AiJ^v that characterize the 
oscillations of E(t) are simply written as 



Em = max{|r|, |rc„t|} 
Ai;7v- = min{|r|, |rcr-it|}. 



(5) 
(6) 



In the above equations r is the initial squeezing factor 
defined in terms of the dispersions of the initial state 
(5x_ and 5p^ as 



r = - In 
2 



5x^ 



and Tcrit is related to the squeezing of the equilibrium 
state for the a;+-oscillator: 



rcrit ^^^^ 
Finally, Scrit is defined as 



A^ 



1 



Scrit = - ln[4Aa;+Ap+(5x_5p_], 



(8) 



(9) 



(7) 



and turns out to be simply related with the entropy of the 
asymptotic state by the symplectic area of the oscillators 
x±. 

Using these results we can conclude that there are three 
qualitatively different types of evolutions for the entan- 
glement for long times. First, entanglement may persists 
for arbitrary long times when ||r| — |rcrit|| > Scrit- In this 
case there is no sudden death of entanglement (NSD). A 
different behavior, an infinite sequence of events of sud- 
den death and sudden revival (SDR), takes place when 
\\r\ - \rcrit\\ < Scrit but \r\ + \rcrit\ > Scru- Finally, a 
third phase characterized by a final event of sudden death 
(SD) of entanglement is realized if |r| + \rcrit\ < Scrit- 



A. Interpretation: Where does the entanglement 
come from? 

The above results may seem, at first sight, somewhat 
puzzling. The final state of the two oscillators may be 
entangled even if there was no entanglement in the initial 
state. In some sense, the common environment provides 
a quantum channel through which entanglement between 
the two oscillators can be either created or destroyed, de- 
pending on the circumstances (initial state, temperature, 
etc). Here, we will discuss this result and present a very 
simple interpretation. The key is to realize that in the 
asymptotic (long time) regime the net effect of the inter- 
action with the environment can be represented by the 
diagram shown in Fig. [l] In the diagram time flows 
from left to right. The original oscillators xi^2 are trans- 
formed into the virtual oscillators x± by the action of 
an ordinary 50/50 beam splitter. After the beam-splitter 
the oscillator x_ evolves unitarily and decoupled from 
which, in turn, interacts with the environment. In 
the long time limit the interaction with the environment 
leads to an equilibrium state for x+ which is completely 
uncorrelated with the state of x_. Such state is Gaus- 
sian and fully characterized by the equilibrium variances 
Aa;_|. and Ap^. Finally, the second beam splitter recom- 
bines the two virtual oscillators to produce again the real 
modes xi and X2- 

The entanglement in the final state is a quantum re- 
source that can be originated from other quantum re- 
sources available in the state of the x± oscillators immedi- 
ately before the second beam splitter. Indeed, such quan- 
tum resource is squeezing. In fact, it is well known that 
a beam splitter can produce entangled states of the out- 
going modes provided the input modes are squeezed ^29"! . 
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5x_,5p_ 




FIG. 1: The evolution of two resonant oscillators coupled to 
a common environment (left side of the diagram) turns out 
to be described, in the long time regime, by the right side of 
diagram: A 50/50 beam splitter combines the original xi^2 
oscillators to form the x± modes. While X- evolves freely, 
x+ couples to the environment which drives this mode to an 
asymptotic equilibrium state which is completely uncorrelated 
with X- and is fully characterized by the dispersions Aa;+ and 
Ap+. A second beam splitter recreates the xi^2 modes which 
will be entangled only if there is squeezing in the x± modes 
before the second beam splitter. 



This observation enables us to understand the origin of 
the entanglement in the final state of the xi^2 oscillators: 
It conies from the squeezing which is available either in 
the X- or in the x+ oscillators. In turn, the squeezing in 
the X- mode, which naturally oscillates with a frequency 
UJ-, is itself inherited from the squeezing (or entangle- 
ment) eventually present in the input modes. On the 
other hand, squeezing in the asymptotic (equilibrium) 
state of the x+ oscillator, which is measured by rcrit, is 
also a source for entanglement in the final state. Below, 
we will show that there are some situations in which a 
non- vanishing value for Vcrit arises as a non-trivial (non- 
Mar kovian) effect of the environment. 

The three phases we mentioned above can be under- 
stood using this interpretation. To make more evident 
the connection between final entanglement and squeezing 
it is convenient to rewrite the equations for the asymp- 
totic entanglement ^ as follows: 

E{t) = \rcrit\ - Scrit + \r\G{t), if \r\ < \rcrit\, 
E{t) = \r\ - Sa„t + \rcr^t\G{t), if \r\ > \r^„t\- 

We can extract some interesting conclusion from these 
equations. First, it is clear that for initial values \r\ < 
\fcrit\ it is possible to use the environment as a resource 
from which we extract entanglement. Thus, for such 
low squeezing the entanglement in the final state can be 
larger than the available quantum resource (squeezing) 
present in the initial state. This is indeed the case if the 
inequality |rcrit| — Scrit ^ 2|r| holds. In other cases the 
environment degrades the quantum resource which is al- 
ready present in the initial state (either in the form of 
squeezing or entanglement). 

With this interpretation, and using the ideas discussed 
in we can conclude that non-classicality at the out- 
put fields (after the second beam splitter) must arise from 
some form of non-classicality at the input. This can ex- 
ist if the equilibrium state has some degree of squeezing 
(which is the case for position coupling) or if the initial 



state of X- is non-classical. For instance, with initial co- 
herent states the condition for the existence of entangle- 
ment in the final state (r = 0, Sx+6p+ = 6x_6p_ = 1/2) 
is rcrit > 1/2 ln(2Ax+Ap+). Thus, to fulfill this con- 
dition we need the environment to produce an equilib- 
rium state where the variance of one of its quadratures is 
smaller than the vacuum limit, i.e. min{A^a:+, A^p+} < 
1/2 (for m= l,f}_ = 1). 



III. MODELS FOR THE ENVIRONMENT 

In this Section we will discuss in some details two 
widely used models for system-environment interaction 
and will discuss the nature of the asymptotic state that 
is obtained. This, as mentioned above, will determine 
the way in which entanglement evolves. We consider a 
system formed by two identical oscillators, with mass m 
and bare frequency luq, whose Hamiltonian is 



^{pI + pD + '^^HxI + xD 



7T!,Ci2Xia;2 -I- Ci2 



~ PlP2 \ 
2 



(10) 



The last two terms include a general type of coupling be- 
tween the oscillators being ci2 and Ci2 the corresponding 
bare coupling constants. The environment will consist 
of a set of harmonic oscillators whose coordinates are la- 
beled as g„. We will consider the coupling between the 
system and the environment as described by 



N 
n=l 



f P4 



N 

E 



^ ^ n— 1 



(11) 



Two cases will be discussed in detail, a) Position cou- 
pling: This corresponds to the case where c„ = 0, i.e. 
the coupling is bilinear in the coordinates a;+ and g„ (in 
this case we will consider that the coupling between the 
oscillators is only through position, i.e. C\2 = 0)- b) Sym- 
metric coupling: This corresponds to the case c„ = c„. In 
this case the interaction is symmetric under interchange 
of position and momentum and can be easily written in 
terms of creation and annihilation operators. It is pre- 
cisely of the form obtained in the so-called RWA (in this 
case we will consider that the interaction between the 
oscillators is also symmetric, i.e. Ci2 = £12). 

For these two cases, the effect of the environ- 
ment is determined by the spectral density defined as 
'^('"^) = '^"^(^ ~ ""^'O/^^nWn- We wiU prcscnt re- 
sults for the so-called Ohmic environment where J(w) — 
^mjQw6{A — w), where A is a high frequency cutoff and 
7o is a coupling constant (other spectral densities were 
analyzed in [2S]). We also assume factorized initial con- 
ditions between system and environment, and that the 
initial state of the environment is thermal (with temper- 
ature T) . For the two types of coupling we consider it is 
possible to obtain an exact master equation that governs 
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the evolution of the reduced density matrix of the two 
oscillators. We will briefly describe them now. 

For position coupling the exact master equation was 
obtained some time ago and reads [3(1: 



- Z?(t)[a;+,[x+,p]]-/(i)[a;+,[p^ 



-p]]- (12) 



Here, the renormalized Hamiltonian is — + 
mSuj'^{t)x'^/2. The effect of the environment shows up in 
four terms: The environment induces a renormalization 
of the frequency of the a;+ oscillator, a dissipative term 
with a time dependent damping rate 7(f) and two dif- 
fusive terms with time dependent coefficients D(t) and 
f{t). All coefficients depend on the environmental spec- 
tral density in a rather complex way {D{t) and f{t) also 
depend on the initial temperature T). The explicit form 
of these coefhcients was studied elsewhere [MtiSl^- Here, 
we will only use the fact that for the Ohmic environment 
all coefficients approach constant asymptotic values in 
the long time limit. It is worth pointing out a tech- 
nical detail concerning the renormalization induced by 
the environment: The bare frequencies of the virtual os- 
cillators are uj^. — ojq ± C12. The interaction with the 
environment renormalizes the frequency of x+ that is 
shifted according to il^{t) = luq + C12 -t- Suj'^{t) (notice- 
ably, the shift Suj'^ diverges in the limit of large cutoff 
A), while the X- oscillator evolves freely with frequency 
= Wg — C12. From these expressions we can also write 
the renormalized frequency and coupling of the real oscil- 
lators n%{t) = Lu^+Suj'^{t)/2 and Ci2(t) = ci2+Su;^{t)/2. 
Therefore for the physical frequencies of both virtual os- 
cillators to be independent of the cutoff one needs to 
absorbe Suj'^ in a renormalization of the bare frequency 
ujQ and the bare coupling constant C12. In such case, 
the long time oscillations of entanglement become cutoff- 
independent. 

The master equation can be used to obtain expressions 
for the variances of position and momentum for the X-^- 
oscillator. In fact, for A^a;+ — {x\) and A^p_|_ = {p\) 
we find 




D 



m ' 



(13) 



and {{x+,p+}) = 0. Where we used upper case letters for 
the renormalized quantities and we omit the time label 
when referring to value of the coefficient in the long time 
regime. 

For symmetric coupling an exact master equation ex- 
ists (for all spectral densities and initial temperatures of 
the environment). In the Appendix we show a simple 
derivation of such equation valid for arbitrary temper- 
atures. Not surprisingly, the equation is nothing but a 
symmetrized version of the previous one (12 1: 



p = -i[HR,p]-i^{t)[[x+,{p+,p}]-[p+,{x+,p}]^ 
- D{t)([x+,[x+,p]] + ^^[p+,[p+,p]]). (14) 



Here, the renormalized Hamiltonian is Hj^ = Hs + 
(5fi^(i)^2^ + ^'^+^+)/'^+- III this case the effect of 
the environment is contained in three terms. The renor- 
malization is also symmetric under position and momen- 
tum interchange. Renormalized mass and frequency of 
the oscillator must be defined as M{t) = m/{l + 

(Sn'^it) + ci2)/lu^) and n{t) = cJo(l + {Sn^{t) + ci2)/^n)- 
As expected, both damping and diffusion terms are sym- 
metric under canonical interchange of position and mo- 
mentum. As in the previous case, all the coefficients of 
the master equation approach constant asymptotic val- 
ues. As before, it is straightforward to obtain the values 
of position and momentum dispersions: 



Ap+ = MnAx+ = 



(15) 



and — 0. One can notice that, contrary to 

what happened with the previous model, the asymptotic 
state of the oscillator x+ has balanced variances. 




FIG. 2: Dynamical phases for the entanglement between two 
non-interacting oscillators, C12 = 0, that are coupled to the 
same environment. Vcrit (black line) and Scrit (grey line). For 
position coupling (left panel) there are three phases (NSD, 
SDR, SD) while for symmetric coupling (right panel) only 
two phases exist. The initial state of the two oscillators is 
such that the virtual oscillator x_ is in a pure state with 
squeezing r. The environment is initially in a thermal state 
with temperature T. 

In |25j the evolution of entanglement was studied for 
these two models under the assumption that the oscilla- 
tors did not interact directly (i.e. C12 = 0) and initial 
conditions such that the state of the x_ oscillator was 
pure (i.e. Sx^Sp^ — 1/2). In such case the values of 
rcrit and Scrit were obtained for a variety of spectral 
densities. The phase diagrams describing the evolution 
of entanglement for the two types of coupling were pre- 
sented and are shown in Figure [2] For position coupling 
the fact that Vcrit may be non-vanishing opens the door 
to the existence of an NSD phase for low temperatures 
and small squeezings. It is also responsible for the exis- 
tence of the SDR phase. These two effects disappear for 
the case of symmetric coupling, that only exhibits two 
phases. Below, we will generalize these results. It is also 
worth pointing out that the above results are valid for res- 
onant oscillators with equal coupling to the environment. 
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for non-resonant oscillators or spatially separeted os- 
cillators |23j . the asymptotic entanglement becomes in- 
dependent of the initial state and sudden death occurs 
above a critical temperature or distance. 



IV. GENERALIZATIONS 

A. Initially mixed states 

As a first generalization we consider the case where 
the initial state of the virtual oscillator x- is mixed. It 
is clear that for this situation all the above formulas ap- 
ply. The only footprint of the initial state of the system 
appears through the dependence of Scrit on the initial 
dispersions of the X- oscillator. The purity of the state 
of X- is characterized by the product 5x-Sp- that ap- 
pears in Scrit- When this product is increased the value 
about which the entanglement oscillates decreases (as can 
be seen from the equations ^ and Q). This implies 
that as a consequence of the impurity present in the x- 
oscillator the final entanglement decreases. This situa- 
tion applies, for instance, to the case in which the initial 
state of each real oscillator is mixed or (using the quan- 
tum optical analogy) for an initial pure state that ends 
entangled after the application of the first beam-splitter. 
The phase diagram for these states can be obtained in a 
simple way from the diagram for initial pure states by a 
translation of the curve Scrit to the right. The net effect 
of this change is to move upwards the horizontal axis. As 
a consequence, the value of the critical temperature be- 
low which the NSD phase exists becomes lower. In fact, 
the NSD island can disappear depending on the degree 
of purity of the state of x^. The phase diagram for both 
types of couplings is shown in Figure [3] 




0.25 0.5 0.75, ,1 1.25 1.5 



FIG. 3: Dynamical phases for the entanglement between two 
non-interacting oscillators initially in a state such that the 
virtual oscillator X- is in a mixed state Sx-5p- = 1, for 
position coupling (left panel) and symmetric coupling (right 
panel). 



of the previous one. When the coupling with the environ- 
ment is through position, the renormalized frequencies of 
the oscillators x± are fij. = Uj^ ± C12, where fi/j is the 
renormalized frequency of the xi^2 oscillators. The non- 
vanishing coupling induces different frequencies for both 
virtual oscillators. In such case the value of Vcrit is: 



rcrit = ^ 1^ 



Ax^ 



1 



In 



IT 



(16) 



The most important difference with the previous case of 
non-interacting oscillators is that rcrit can be different 
from zero even in the limit where the x+ oscillator is not 
squeezed. As a consequence, it is possible to observe os- 
cillations of the entanglement in the high temperature 
regime. Some examples of the behavior of the entan- 
glement in this situation are shown in Fig. |4] (a). It 
is interesting to stress that, as we mentioned above, in 
all cases it is necessary to include a bare coupling term 
proportional to Suj'^ in the Hamiltonian. Thus, only in 
this way the system would not oscillate with a cutoff- 
dependent frequency in the long time limit (see Fig. |4] 
(6) and (c)). 




time (a. u.) 



FIG. 4: Dynamics of the entanglement for resonant oscillators 
and position coupling at T/fijj — 10, 70 = 0.1 and A — 20. 
(a) Initial two-mode squeezed state (r — 3, f2i,2 = i^R = 1) 
with different couplings, we observe oscillations at high tem- 
peratures when there is non-vanishing renormalized coupling 
C12. (6) Initial two-mode squeezed state (r = 3, — 3), 
oscillations that depend on the cutoff frequency are present 
when one consider a vanishing bare coupling (i.e. C12 — and 
C12 ~ 5lo'^/2). (c) Separable initial state (r = 3), C12 ~ 0. 
The final entanglement depends on the degree of purity of the 
initial state and is lower than the entanglement achieved for 
an initial pure state in a quantity given by In [2Sx-Sp-] /2. 



B. Interacting oscillators 

If the renormalized coupling C12 = C12 -I- Suj'^/2 is non- 
vanishing the analysis is also a straightforward extension 



For the case of symmetric coupling the asymptotic be- 
havior does not change considerably if we add an inter- 
action between the oscillators (provided the interaction 
is also symmetric in position and momentum). In fact, 
in this case we always obtain Mfi = m-Uj_ = muQ and 
a vanishing rcrit , as it can be seen from the equations (l8| 
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and (15 1. As a consequence, for the symmetric coupling, 



entanglement is constant in the asymptotic regime and 
the phase diagram does not change. The phase diagram 
for position coupling with interacting oscillators is shown 
in Figure [5j 



T/Q. 




FIG. 5: Dynamical phases for the entanglement between 
two interacting oscillators, with C12 = —0.5, that are cou- 
pled to the same environment. The initial state is such that 
Sx-Sp- = 1/2. Vcrit achieves a constant value at high tem- 
peratures, and a SDR region appears at high temperatures. 



V. CONCLUSIONS 

In summary, we have presented a simple general 
overview of the evolution for the entanglement between 
two oscillators coupled to a common Ohmic reservoir. 
We showed how the existence of non-trivial phases for 
the evolution of the entanglement can be seen by us- 
ing a simple interpretation based on a quantum optical 
analogy. The fact that the x± virtual oscillators decou- 
ple and that the x+ oscillator approaches equilibrium is 
what makes this interpretation possible. The nature of 
the equilibrium state for a;+ may be peculiar since it may 
be squeezed by a factor Vcrit when the coupling with the 
environment is not symmetric under position-momentum 
interchange. In the final Section, we generalized the 
results previously obtained in |M1 US] to include cases 
where the state of the virtual oscillators is mixed as well 
as to consider the case of interacting systems. In all cases, 
the evolution of entanglement can be discussed in terms 
of a phase diagram that describes all qualitatively differ- 
ent behaviors in the long time limit. It is worth noticing 
that the phase diagram for interacting oscillators shown 
in Figure |5] that includes the three dynamical phases 
(NSD, SDR and SD), seems to be observable in exper- 
iments realizable with current technologies in ion traps 



Acknowledgments 

JPP is a member of CONICET and AR has a fel- 
lowship from CONICET. This work was supported with 
grants from ANPCyT (Argentina) and Santa Fe Institute 
(SFI, USA). 



APPENDIX A: EXACT MASTER EQUATION 
FOR SYMMETRIC COUPLING IN POSITION 
AND MOMENTUM 

We present here a brief sketch of a simple derivation of 
the exact master equation for symmetric coupling which 
is similar to the one obtained for position coupling in 
[53 . The derivation is valid for all spectral densities and 
temperatures. The full Hamiltonian, written in terms of 
creation and annihilation operators is: 



H = huja^a 



+ ynk{abl + bal). (Al) 



The first step of the derivation is to notice that this 
Hamiltonian preserves the Gaussian nature of the states 
(this is obvious since is a quadratic form of the co- 
ordinates and momenta of all the particles). Therefore, 
the evolution operator of the reduced density matrix of 
the system is a Gaussian operator also. It is possible to 
show, following the same steps described in the deriva- 
tion contained in [SI] , that if the propagator is Gaussian 
the master equation is local in time. Moreover, it has a 
limited number of terms whose number is further reduced 
if one imposes the constraint that the equation should be 
symmetric under canonical exchange of position and mo- 
mentum. Thus, under this condition one can show that 
the form of the master equation should be 



+ 7(f) (2apa^ — a^ap — pa} a) 



"f(t) (2a) pa — aa^ p — paa^) 



(A2) 



This equation contains three unknown coefficients with a 
clear physical interpretation: A renormalized frequency 
in Hi^{t) such that rifl,(t) = lu + (5f2^(i)/w^, a dissipa- 
tion coefficient j{t) and a diffusion coefficient D{t). The 
above argument simply tells us that the master equation 
should have this form but does not enforce any constraint 
in the time dependence of such coefficients. Now, we will 
find them using the following argument. 

From the total Hamiltonian, we can derive Heisenberg 
equations for all the operators, which turn out to be: 



da 
lit 



k 



(A3) 
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dbk 
dt 

dt 



= iwkhk - igka, 



These can be formally solved as 

a{t) ^ ua{0) + Y,Pnbn{0), 
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bk{t) = dka{t) + y^gfcn5„(0), 



(A4) 
(A5) 

(A6) 
(A7) 



where u, p„, dk and qkn are appropriate time-dependent 
coefficients. 



On the other hand, the master equation (A2| can be 



used to obtain evolution equation for expectation values 
of the operators of the system which turn out to be: 



d{a) 
"dT 



— , = -A^{t){aa) + a) a) + 

dt muj 



(A8) 
(A9) 



Comparing equations (A3 1 with the expectation value 



of equation ( A5 ) we can simply obtain the expressions 



for the time dependent coefficients. Moreover, impos- 
ing that the initial state of the environment is thermal, 
{b\{0)bk{0)) = Uk, these coefficients can be expressed as: 



7(t) 

sn^t) 

m 



^^9k{dk- dl), 

k 

^^9k{dk + dl), 

k 

j^9k{q*kiPi - qkiP*i){2ni + 1). 



(AlO) 



Hence, they are completely defined in terms of the solu- 



tion to the equation of motion (A6l-(A7l 



The solutions a bove equations can be explicitly writ- 
ten. In fact, from (A4| we can write: 



bk{t) = e-™'=*6fe(0) - igk / e-'^'^'-'^ais)ds. (All) 



Moreover using (A3 1 and (All I we get 



da{s) 



+ iuja(s) + / ri(s — x)a{x)dx = if(s), (A12) 
as Jq 



where /(s) = — J^k 9^^ ™'°*^/c(0); and the kernel r/(s) is 
defined as 



77(5) = / dwJ{w)e 
Jo 



(A13) 



Here, the spectral density is J{w) — X]fe5fc'^(^ ~ Wk)- 



Now, the equation (A12) can be solved using the Laplace 
transform and the convolution theorem, with initial and 
final conditions given by: a{s = 0) = a(0) and a(s = t) = 
a{t). From these expressions it would be possible to find 
the time-dependence of all the coefficients of equations 
(A6l-(A7l. Finally, we note that there are certain rela- 



tions satisfied by the coefficients: One of them arises from 
the fact that the commutation relations are preserved 
during the evolution. Thus, as [a{t),bk{t)^ = 0, one 
has dk = — J2n QknPn- This relation can be used to show 
that at zero temperature we have D{t)/mLjj = 7(i). And 
the exact master equation has only two time-dependent 
coefficients: frequency renormalization and dissipation. 



k.l 
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